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The dynamics of many glassy systems are known to exhibit string-like hopping motions each 
consisting of a line of particles displacing one and other. By using molecular dynamics simulations 
of glassy polymers, we show that these motions become highly repetitive back-and-forth motions 
as temperature decreases and do not necessarily contribute to net displacements. Particle hops 
which constitute string-like motions are reversed with a high probability, reaching 73% and beyond 
at low temperature. Structural relaxation rate is then dictated not by a simple particle hopping 
rate but instead by the rate at which particles break away from hopping repetitions. We propose 
that disruption of string repetitions and hence also structural relaxations are brought about by 
pair-interactions between strings. 


I. INTRODUCTION 

When cooled sufficiently rapidly, a liquid can be 
quenched into the glassy state, where motions of the 
constituent particles become dramatically slowed down 
resulting in solid-like responses. Despite its prevalence, 
glass transition remains poorly understood m- Fun¬ 
damental questions including the nature of the transi¬ 
tion and the mechanism for the dramatic slowing down 
of the dynamics are still not settled. Leading theories of 
glassy dynamics include Adam-Gibbs theory [4], random 
first order transition (RFOT) theory [5j [6], elastic theory 
|7l, mode-coupling theory [8], dynamic facilitation jMT2] , 
and so on. 

Molecular dynamics (MD) simulations allow access to 
the microscopic details of the dynamics m • A remark¬ 
able discovery is string-like motion in which neighboring 
particles arranged in a line displacing one another from 
head to tail m- They have been observed in simula¬ 
tions of fragile glasses including binary Lennard-Jones 
liquids 

PUCES], monoatomic Dzugutov liquids m , and 
bead-spring polymers m- They are also observed in 
experiments of colloidal systems T5J [19] . String-like mo¬ 
tions in general are incoherent and consist of independent 
micro-strings each of which makes simultaneous coherent 
hops [16]. Micro-strings are more bounded in length and 
most are of 1 to 3 particles long m ■ String-like motions 
are often considered as structural relaxations, i.e. a re¬ 
laxations m • Micro-strings have been observed to be 
reversible [T6j 20 and have thus been considered as (3 re¬ 
laxations m ■ String-like motions extracted from short- 
time displacements, which are not necessarily coherent, 
have also been interpreted as (3 relaxations [22]. How¬ 
ever, detailed criteria and mechanisms concerning how (3 
relaxations are related to structural relaxations is lack¬ 
ing. 

String-like motions form the basis of several interest¬ 
ing theories of glassy dynamics. Most of them assume 
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the existence of a term in the structural relaxation en¬ 
ergy barrier proportional to the average string length l 
[6j [23j [24]. Assuming also a significant increase in l as 
temperature T decreases as suggested in simulations im 
kinetic arrest follows. 

An important recent advancement by Chandler and 
coworkers EUHg is the suggestion of a physical realiza¬ 
tion of hierarchically constrained dynamics m based on 
string-like motions. The authors identified an elementary 
excitation as a structural feature allowing a micro-string 
motion and a larger particle displacement indicates an 
excitation at a higher level. A simple dependence of the 
dynamic rate on level was then suggested. The detailed 
mechanism dictating how each hierarchical level facili¬ 
tates or constrains each other deserves further investiga¬ 
tion. 

Glasses are routinely studied in both equilibrium and 
non-equilibrium states m mi- Besides investigating 
equilibrium samples of a glassy polymer system from 
MD simulations, we also examine non-equilibrium deeply 
supercooled samples quenched to much lower T. We 
show that as T decreases, string-like particle hopping mo¬ 
tions dominate the dynamics and hopping distances be¬ 
come concentrated around the particle diameter. We find 
that hoppings and thus also string-like motions become 
increasingly reversible and repetitive. We argue that 
string-like motions, as opposed to only micro-strings, are 
mostly (3 relaxations at low T. Since identical string-like 
motions can be reversed and repeated multiple times, we 
refer to them as multiple back-and-forth propagation of 
a single string. We suggest that the main mechanism 
rendering string-like motions less readily reversed and 
may hence contribute to structure relaxations is pair- 
interactions between strings. Pair-interactions also de¬ 
stroy old strings and facilitate the creation of new strings, 
resulting in effective string mobility. Examples of string 
repetitions and interactions are illustrated by particle 
trajectories from MD simulations. Systematic quanti¬ 
tative assessment of string repetition based on particle 
hopping statistics are performed. 

The rest of this paper is organized as follows. Sec¬ 
tion in describes our simulation method. Section mu- 
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lustrates string repetition. Quantitative measurements 
on particle hoppings and their repetitive properties are 
presented respectively in Sec. IV and [V] We explain in 
Sec. |VI| the disruption of string repetitions by string 
interactions. Section |VII| discusses the relationship be¬ 
tween voids, strings and mobile clusters while Sec. \vm\ 


concludes the paper with a summary and some further 
discussions. 


II. MODEL AND SIMULATION METHOD 


Our simulations are based on a bead-spring model of 
polymer following Refs. I28H30] generalized to heavier 
chain-tails for producing a uniform monomer mobility. 
Specifically, pairs of particles interact via the Lennard 

Jones (LJ) potential 4e (cr/r) 12 — ( a/r ) 6 with a cutoff 

2-2 1 / 6 beyond which it becomes a constant. The potential 
hence has both attractive and repulsive parts. Bonded 
particles are further bounded by a finitely extensible non¬ 
linear elastic (FENE) potential — |Rq In 1 — (r/Ro) 2 

where k = 30 e/a 2 and Rq = 1.5cr. We adopt dimension¬ 
less LJ units which amounts to taking a = e = 1. We 
nevertheless express lengths in unit of a for clarity. Inter¬ 
nal particles in each chain have a mass M/ = 1. However, 
in contrast to common practice [28H30] . one particle at 
each end of the chain has a larger mass Mt = 4. This 
leads to a uniform mobility for all particles in the chain 
as is verified from measured particle mean square dis¬ 
placements (MSD). 

Our main simulations involve two samples of polymer 
melts each having 1500 chains of length 10 inside a cubic 
box with periodic boundary conditions simulated using 
a timestep of 0.005. The melts are first thermalized at 
T = 0.6 and then 0.5 following standard techniques [28]. 
They are then repeatedly quenched and annealed under 
NPT conditions at zero pressure by steps of AT = 0.02 to 
produce samples at various T. Each of these quenching 
or annealing processes involves 10 7 timesteps. A sample 
is further annealed for 10 8 timesteps before data taking. 

We report our main studies at temperature T in be¬ 
tween 0.30 to 0.40 exhibiting dynamics ranging from 
those of deeply to slightly supercooled liquids. A sample 
is considered equilibrium if the particle MSD during the 
final annealing process reaches 5cr 2 . For T = 0.36, 0.38 
and 0.40, fully equilibrated polymer samples are stud¬ 
ied. For T = 0.30, 0.32, 0.34 and 0.36, we consider non- 
equilibrium samples quenched from higher T. Note that 
at T = 0.36, we study both equilibrium and quenched 
samples for comparison. Specifically, the sample prepa¬ 
ration procedure described above leads to the equilibrium 
samples for T = 0.38 and 0.40 and the quenched samples 
at T = 0.30,..., 0.36. For T = 0.36, we further perform 
an extended annealing of 10 10 timesteps and this gener¬ 
ates two additional samples at T = 0.36 satisfying the 
above equilibrium criterion. 

All our analysis are based on coarse-grained particle 


trajectories r- c (£) recorded during data collection runs 
each consisting of 10 9 additional simulation timesteps un¬ 
der NVT conditions. Specifically, each value of r^ c (£) is 
a coarse-grained position of particle i defined by 

)t / e[t,t+Atj (1) 

averaged over a duration A t c = 5 (i.e. 1000 timesteps). 
In practice, this is evaluated by an in-situ average over 
20 instantaneous positions fi(t') of particle i taken af¬ 
ter every 50 timesteps. The coarsening time A t c = 5 is 
long enough to average out most vibrations about meta¬ 
stable positions. It is also negligible compared with the 
typical waiting time between two consecutive hops of 
a particle. Therefore, r^ c (£) nearly always points to a 
meta-stable particle position, rather than somewhere in¬ 
terpolating between two meta-stable positions related by 
a hop. This improves the accuracy of our measured hop¬ 
ping statistics. For each sample, we record 10 4 snapshots 
of coarse-grained positions r^it) taken after every dura¬ 
tion of St = 500 (i.e. 10 5 timesteps). Note that St is still 
much shorter than the average waiting time between par¬ 
ticle hops so that motions neglected in between snapshots 
only have small quantitative impacts on our results. 

Quantitative measurements to be presented are av¬ 
eraged over two independent samples, for which error 
bars are generally smaller than the symbol sizes in our 
plots. They are based on r^ c (£) and are in general consis¬ 
tent with those using fi(t) except at length scales much 
smaller than a where the reduced fluctuations may lead 
to quantitative impacts. Simulations are executed us¬ 
ing the highly efficient HOOMD package m on Nvidia 
GTX580 GPU’s. Each extended annealing at T = 0.36 
takes 5 months of runtime. Our MD simulations at 
T = 0.36 each involving more than 10 10 timesteps is in 
our knowledge far more intensive than similar previous 
studies. Equilibration at T < 0.34 will be substantially 
more challenging computationally and should be imprac¬ 
tical with current computing technology. 


III. STRING REPETITION 

At low temperature T, particle motions mainly consist 
of thermal vibrations and occasional hops. Figure [lja)- 
(b) shows a typical hop of particle i from a quenched 
sample at T = 0.32 as represented by the true trajectory 
fi(t) recorded at every MD timestep [32] and the coarse¬ 
grained trajectory rfit) explained in Sec. |TI[ As seen, the 
coarse-grained trajectory clearly illustrates the hopping 
motion and will be adopted throughout this work. We 
consider particle i hopped during a given period if the 
recorded r^ c (£) shows any displacement beyond a distance 
of 0.6cr within any sub-interval of the period. 

It is instructive to visually identify string-like motions 
directly by plotting the trajectories of all the hopped par¬ 
ticles in a neighborhood. Figure [ljc) shows a typical 
string-like motion. Note that the particles involved in 
general do not belong to the same polymer chain m- 
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FIG. 1. Particle trajectories colored according to time t with 
red (blue) being the earliest (latest), (a) Hopping of a single 
particle as shown by its true trajectory ri(t). The particle 
first vibrates about one site and then hops to another site to 
the right, (b) The same motion as in (a) represented by the 
coarse-grained trajectory consisting of line segments joining 
time-averaged positions f^ c (t). Vibrations before (red) and 
after (blue) the hop are deemphasized by the averaging. The 
long line segment (green) of length about 0.9cr implies a sud¬ 
den large displacement, i.e. a hop. (c) Coarse-grained trajec¬ 
tories of 3 particles arranged in a line performing a string-like 
motion. Particles hop to the right and nearly fill the positions 
of the preceding particles. The uniform color (green) of the 3 
long segments implies simultaneous hops, (d) Coarse-grained 
trajectories of the same particles in (c) shown for a longer 
duration. Four string-like motions (red, orange, light blue, 
blue) occur, the first of which is the one in (c). The durations 
shown are 200 in (a)-(b), 5,000 in (c) and 20,000 in (d), while 
T = 0.32 in LJ units. 


A key observation that has motivated this work is that 
string-like motions become highly repetitive exhibiting 
strong memory effects at low T. For the string-like mo¬ 
tion shown in Fig. ljc), it is subsequently reversed, re¬ 
peated and reversed again as shown in Fig. 0d)- There 
are thus 4 successive string-like motions. 

Figures [2] and [3] show time-consecutive coarse-grained 
trajectories of all hopped particles from an equilibrium 
sample at T = 0.36 and a quenched sample at T = 0.32 
respectively. Line segments constituting the trajectories 
are colored according to their positions in lateral direc¬ 
tion. This allows an indication of depth, as spatially 
overlapping trajectories of distinct colors must be well 
separated in three dimensions. We can observe dozens 
of examples of string repetition. We have also produced 
movies showing the evolution of the trajectories for much 
longer durations [33] . 

From various visualizations, we observe repetitions of 
both coherent and incoherent string-like motions of dif- 




FIG. 2. Coarse-grained trajectories of hopped particles dur¬ 
ing three consecutive durations each of 2500 from an equilib¬ 
rium sample at T — 0.36. Displacements are colored from red 
to blue starting from the left surface of the cubic simulation 
box of width 24.3a. 


ferent lengths. Incoherent strings can be either fully or 
partially reversed and repeated. We also provide an ani¬ 
mation of the repetition of an incoherent string-like mo¬ 
tion [33] . Persistent repetitions over long durations usu¬ 
ally occur to those far away from the other strings. Rep¬ 
etitions become increasingly common and robust as T 
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FIG. 3. Coarse-grained trajectories of hopped particles dur¬ 
ing three consecutive durations each of 5000 from a quenched 
sample at T = 0.32. The coloring is similar to that in Fig. [2] 

decreases. 

String-like motions at higher T are often considered 
structural relaxations m- Observing here that they 
may not contribute to net displacements, these oscilla¬ 
tory parts of the motions should be classified only as 
caging dynamics or /? relaxations alongside thermal vi¬ 
brations. We will explain later that the dominant mech¬ 
anism to break away from the indefinite repetitions of 


string-like motions is a pair interaction between strings. 
After repetitions finally stop and assume an odd number 
of hops, only the last outstanding hop of each particle 
contributes to structural relaxations. 

String-like motions can be easily identified visually 
from particle trajectories obtained from MD simulations 
as shown above or from potential energy landscape (PEL) 
calculations [34]. Conventionally, they are often ex¬ 
tracted using standard methods which also dictate their 
precise working definitions El- Algorithmically, the net 
displacements of particles over a time period of At com¬ 
parable to an average structural relaxation time are typ¬ 
ically studied. Particles with the top 5% displacements 
during At are grouped into clusters of adjacent mobile 
particles. Each cluster is then resolved into string-like 
motions on the basis that particles belonging to a string 
best displace the preceding members. 

We have also identified and quantitatively analyzed 
string-like motions following conventional algorithms PH 
adapted to considering particle displacements beyond 
0.6cr over relatively short time periods At to allow appli¬ 
cations to low T. String lengths are found to nicely follow 
an exponential distribution down to the string length of 
1 particle with an average string length l ~ 2 particles 
varying only mildly with At and T. This is consistent 
with recent results on Ni nano-particles [22] . By defi¬ 
nition, all hops belong to some string-like motions as a 
string can be one particle long. In the following, we will 
simply focus on particle hopping statistics which have 
straightforward implications on string statistics. 

IV. HOPPING DISTANCE 

We define particle displacement r over time r based on 
coarse-grained trajectories as 

r = | G c (t 0 + t) - G c (£ 0 ) | . (2) 

We calculate the probability density P T (r) of r, which is 
often written as P T (r) = 47rr 2 G s (r, r) using the van Hove 
self-correlation function [26] 

1 N 

G a (r, t) m 5>( W + t) - - r)) (3) 

V i= 1 

where N is the total number of particles. Figure 
plots 47rr 2 G s (r, t) for T = 0.40. At small r, it shows 
a main peak at r ~ O.lcr. It corresponds to small dis¬ 
placements resulting from both structural relaxations as 
well as vibrations that are not completely eliminated by 
the coarsening process. As r increases, the peak broad¬ 
ens gradually as expected from relaxations induced by 
liquid-like collective flow. The distribution however de¬ 
viates significantly from a Gaussian with a weaker sec¬ 
ondary peak at r ~ 0.9a resulting from activated particle 
hopping like those shown in Fig. [l] Similar secondary 
peaks have long been observed for binary Lennard Jones 
liquids [HI HUES and glassy colloidal liquids [35] 
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FIG. 4. Semi-log plot of probability density P T (r) — 
47t r 2 G s (r, t) of particle displacement r during time r at tem¬ 
perature T — 0.40 (a) and 0.36 (b) measured from equilibrium 
samples. 


Figure Qb) shows 47 t r 2 G s (r, t) at T = 0.36 for the 
fully equilibrated samples. The secondary peak at 0.9a 
is now much more prominent. As r increases, it gains in 
height dramatically. In contrast, the main peak broadens 
relatively mildly over the same range of r. The impor¬ 
tance of hopping as compared to collective flow is thus 
significantly increased at this lower T. Note that hopping 
also induces broadening of the main peak as neighboring 
non-hopping particles are displaced due to longer-ranged 
elastic distortions. Therefore, collective flow may have 
diminished even more than naively suggested by the slow 
broadening of the main peak. 

Figure [H] compares 47rr 2 G s (r, r) at various T with r 
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FIG. 5. Semi-log plot of probability density P T {r ) = 

4tv r 2 G s (r, r) of particle displacement r during time r cor¬ 
responding to a MSD of 0.09cr 2 . The dashed line indicates a 
fit to the main peak for the equilibrium case at T — 0.36. 


chosen so that the MSD in all cases are standardized 
at 0.09a 2 , a value attainable for our data at T > 0.32. 
From the equilibrium samples, we observe a monotonic 
strengthening of the secondary peak as T decreases im¬ 
plying enhanced importance of hopping as compared 
to collective flow. Comparing between equilibrium and 
quenched cases at T = 0.36, we find that equilibration en¬ 
hances the relative importance of hopping. The quenched 
samples at lower T however show a weaker T dependence 
and this may indicate that collective flow is close to being 
completely suppressed. 

We further study the equilibrium case at T = 0.36 
in detail. The main peak can be fitted to the empir¬ 
ical form f(r) = 200r 2 (260e _50r + 200 -15r ) shown by 
the dashed line in Figure [5j This portion of the distri¬ 
bution contributes f r 2 fdr ~ 0.0103cr 2 to the MSD and 
accounts for collective flow, elastic distortion due to par¬ 
ticle hops, and vibrations not completely coarsened out. 
The complementary distribution 47rr 2 G s — / capturing 
the secondary peak in contrast contributes 0.0797a 2 to 
the MSD and accounts for particle hops. As a result, 
particle hops are responsible for at least 87% of the MSD 
and clearly dominate the dynamics. This dominance in 
the MSD is even stronger at larger r as can be more read¬ 
ily observed from plots of 47rr 4 G s (r, r) [20] . We therefore 
believe that particle hops of distances centered around a 
completely dominate structural relaxations at low T. 

From Fig. [4]and Fig. [5j the dominant hopping distance 
of about 0.9a is quite independent of r and T. Being 
close to a, it is also close to the distance between meta¬ 
stable particle positions and is thus favored. A further 
subtle point is that being slightly smaller than a, a par- 
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tide usually slightly fall short of reaching a neighboring 
metastable site as revealed by a gap between trajecto¬ 
ries of adjacent particles in string-like motions (see Fig. 
[ljc)). This may indicate some elastic restoring forces 
counteracting the hop. Finally, 0.6cr marks the approxi¬ 
mate position of the dip which is also quite independent 
of r and T. It is adopted as the threshold displacement 
for defining hops as used above and throughout this work. 
In contrast, we find that the conventional approach of 
inferring a threshold from the top 5% most mobile par¬ 
ticles na may not always best distinguish between the 
two peaks observed here. 

V. HOPPING RATES AND RETURN 
PROBABILITY 

We calculate the particle hopping rate R\ = Pi/St 
using the probability Pi that particle i hops during 
[to, to + St] where St = 500. Figure [6](a) plots Ri against 
1/T. The equilibrium data is consistent with a slightly 
super-Arrhenius behavior. The non-equilibrium data 
show more moderate decrease of Ri as T is decreased 
further. 

We now focus on the subsequent motion of a hopped 
particle i by evaluating the probabilities P ret and P 2 that 
it will first perform, respectively, a returning or non¬ 
returning second hop. Specifically, after the first hop dur¬ 
ing [to, to + St ], a second hop is defined as returning if the 
particle returns to within a distance of 0.3cr from the orig¬ 
inal position r^ c (to). A non-returning second hop means 
that it has displaced again elsewhere beyond a distance 
of 0.6cr from the hopped position r- c (to + St). We moni¬ 
tor the particle over a long duration of 10 6 after the first 
hop so that subsequent motion can be categorized in most 
cases. Figure ^b) plots P re t and P 2 against 1/T. As T 
decreases, P ret rises while P 2 drops implying increasing 
memory effects. We find that P ret + P 2 > 0.95 so that 
uncategorized cases are few. These non-returning second 
hops are responsible for the broader peak of 47rr 2 G s (r, r) 
at r ~ 1.6cr as observable from Fig. [5] 

More importantly, P ret reaches a high value of 0.73 
from our equilibrium data at T = 0.36. It is well beyond 
the value of 0.21 found at T = 0.40 more akin to simple 
diffusion. For the quenched samples at T = 0.36, P ret 
takes a slightly smaller value of 0.65, which is expected 
as quenched systems are more driven to previously un¬ 
visited states in attaining equilibrium. For T < 0.34, 
the non-equilibrium values of P ret shown in Fig [6](b) are 
hence likely the lower bounds of the equilibrium ones 
which cannot be directly measured. We find that P ret 
continues to increase as T decreases and reaches a very 
high probability of 0.93 when we extend our study to 
T = 0.26 (data not shown). This leads us to postulate 
that for both equilibrium and quenched samples, 

Pret —^ 1 and P 2 —y 0 (4) 



1/T 



f/T 

FIG. 6. (a) Plot of particle hopping rate Pi, rate of non¬ 
returning second hop P 2 , and MSD-based structural relax¬ 
ation rate R a vs 1/T for T = 0.40,0.38,..., 0.30. (b) Plot 
of probabilities of returning hop P re t and non-returning sec¬ 
ond hop P 2 . In (a) and (b), results based on equilibrium 
(quenched) samples are represented by solid (open) symbols. 

Figure^ a) also plots the rate P 2 of non-returning sec¬ 
ond hops defined by 

R 2 = RiP2- (5) 

The equilibrium data shows that P 2 is slightly super- 
Arrhenius, like Pi. However, due to the hopping repe¬ 
titions, P 2 clearly drops faster than Pi as T decreases. 


as T —0. 
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FIG. 7. Displacement pathways of 6 particles during a pair- 
interaction between two strings illustrated using (a) coarse¬ 
grained trajectories from a MD simulation and (b) the cor¬ 
responding schematic diagram. Trajectories in (a) covering 
a duration of 4500 are colored according to time similar to 
those in Fig. [l] The two distinct colors (green and blue) of 
the 7 long segments imply two set of hops performed at dif¬ 
ferent times. The first string (green) involves the hopping of 
3 particles. It pushes particle j to a new position along the 
second string (blue) whereby the 4 particles the second string 
contains can hop subsequently. Notice that particle j, which 
is shared by both strings, has hopped twice while each of the 
other 5 particles have hopped only once. The temperature is 
0.32. 


This trend is also demonstrated by the non-equilibrium 
data at lower T. 

In addition, we have measured a MSD-based structural 
relaxation time r a [36] which is the time it takes for the 
particle MSD to reach (0.6cr) 2 . This leads to a structural 
relaxation rate defined by 

Ra = 1 /r a . (6) 

We plot R a in Fig. |6|a) for T > 0.34 at which the MSD 
is able to reach (0.6o) 2 during the data collection runs. 

For T > 0.38, Fig. [6ja) reveals f?i ^ f? 2 - R a im¬ 
plying that most hops constitute structural relaxations. 
In sharp contrast, for T < 0.36 for both equilibrium 
and quenched cases, Ri R a . It shows that particle 
hops, which have become increasingly repetitive (as in¬ 
dicated by P ret in Fig. [6^b)), are more frequent than 
the rate of structural relaxations. This supports the as¬ 
sertion that oscillatory hops are (3 relaxations. On the 
other hand, R 2 — R a still holds. Therefore, the second 
non-returning hops, which break string repetitions, now 
dictate the structural relaxation rate. 


VI. STRING INTERACTION 

To understand the mechanism responsible for the non¬ 
returning second hops, we have visually examined hun¬ 
dreds of relevant particle trajectories. Interpretation for 
their origin get significantly easier at low T as string¬ 
like hopping motions completely dominate over collective 
flow. For T < 0.32, all non-returning second hops can be 
identified as being caused by pair-interactions of strings 
similar to the example of string crossing shown in Fig. 


0a) -(b). Here, a particle labeled as j is shared by two 
strings. The first string-like motion (green) induces the 
first hop of particle j. The new position of j enables a 
new string (blue). The second string-like motion along 
the new string then induces the second hop of particle j. 
This however disable the first string so that it cannot be 
reversed unless the second string is first reversed. If not 
reversed eventually, all 7 hops contribute to structural 
relaxation. 

We thus conclude from our observations hat string-like 
repetitions at low T occur indefinitely exhibiting memory 
effects and contribute only to /3 relaxations. Pair inter¬ 
actions can break the repetition and generate structural 
relaxations. At slightly higher T, this mechanism can 
still be observed but are sometimes complicated by the 
presence of liquid-like collective flow. 

Note that conventional algorithms m outlined in Sec. 
Ill fail to properly resolve the strings in Fig. [7|^a). Con¬ 


sidering only the net displacement, conventional algo¬ 
rithms overlook the turn made by particle j. With the 
turn straightened up, particle j is instead assigned ac¬ 
cording to the best-displacing rule to an upward-going 
string of 4 particles on the right. The other 2 parti¬ 
cles constitute another upward going string on the left. 
The algorithms thus report two non-intersecting incoher¬ 
ent strings. This hides the detailed time sequence of 
the dynamics and the phenomenon of string interaction. 
We have found it straightforward to identify any specific 
string in ways consistent with the dynamics by visually 
inspecting time-colored particle trajectories taken over 
some manually selected time-windows, but a practical 
algorithm remains to be worked out. 


VII. VOID, STRING, CLUSTER AND CLUSTER 
REPETITION 

Our understanding of the relationship between conven¬ 
tionally defined strings and mobile clusters m based on 
close examinations of particle trajectories from MD sim¬ 
ulations is as follows. For a specific period At, a number 
of string-like motions and pair interactions may have oc¬ 
curred. Using conventional algorithms that examine only 
the net displacements of hopped particles, strings prop¬ 
agated once or repeated an odd number of times lead 
to the observations of the associated string-like motions. 
On the other hand, particles in strings repeated an even 
number of times are classified as immobile since they re¬ 
turn to the original positions. The associated string-like 
motions are not observed. Furthermore, pair interac¬ 
tions of strings can lead to cross-links between string-like 
motions. Those observed string-like motions in general 
cross-linked constitute the mobile clusters. 

As explained above, string-like motions become in¬ 
creasingly repetitive as T decreases and repetitions can 
be broken by string interactions. In fact, as T further 
decreases, we observe that pair interactions can also be 
reversed and repeated. More generally, we observe rever- 



sal and repetition of sequences of string-like motions and 
pair-interactions, and thus also of mobile clusters. An 
example of a partially repeated cluster of strings shaded 
in blue can be seen at the lower right corner of the sim¬ 
ulation box in Fig. [3] Here, we note that the repetitive 
nature of strings and more generally clusters are due to 
the momentarily frozen self-generated disorders charac¬ 
teristic of glassy fluids. The recurring clusters exemplify 
memory effects at longer time scales. Structural relax¬ 
ations then require breaking the repetition of clusters by 
larger collections of string interactions. 

Free volumes are long known to be important in glassy 
dynamics m- Due to hard-core intermolecular poten¬ 
tials, local spikes in the particle density are unlikely and 
free volumes characterize most density fluctuations. A 
contiguous free volume comparable to the average parti¬ 
cle volume H can be called a void. A particle hop of a 
distance about a is equivalent to the hop of a void over 
the same distance in the reversed direction. However, it is 
easy to see from simple geometry that a shorter hop of a 
particle does not mean a shorter hop of a void. Instead, a 
particle hop of displacement say 0.5a corresponds to the 
hop of half a void over a distance a. Since hops are domi¬ 
nated by distances of about a as shown in Sec. [TV| voids 
are transported mostly in whole. Voids thus provide a 
useful decription of the dynamics we have observed. 

Nevertheless, static particle configurations show that 
voids of volume U are rare. Smaller voids each of vol¬ 
ume a fraction of U are observed in simulations but are 
correlated only slightly to hopping energy barriers [ 54] . 
Voronoi volumes of individual particles have been used 
to indicate free volumes but they also show only small 
correlations with motions [381440] . 

To understand the apparent conflict about the rele¬ 
vance of voids as suggested by the dynamics versus the 
statics, we have examined particle trajectories around 
string-like motions in detail. Since we define particle 
hops with a threshold displacement 0.6cr, strings like that 
in Fig. [ljc) terminates abruptly. More detailed exami¬ 
nations indicate that on each end of a string, there are 
typically a few neighboring particles with smaller dis¬ 
placements decreasing continuously further away from 
the string end. The free-volume H transported by a hop is 
hence distributed across a few interstitial sites before and 
after the hop. From static observations, their existence 
may be easily overshadowed by structural disorder. De¬ 
tection using Voronoi particle volumes can be even more 
challenging as each interstitial volume is further shared 
among the surrounding particles. In addition, a free vol¬ 
ume may induce a hop and then be transported away, 
further weakening the spatial correlation between parti¬ 
cle hops and free-volumes. 

Therefore, we define a free volume H transported by 
a hop and distributed among a few neighboring inter¬ 
stitial sites as a fragmented void. It is a spatially non¬ 
contiguous quasi-particle of a combined volume H dom¬ 
inated by a small void a fraction of H at one end of a 
micro-string accompanied by a few even smaller voids 


nearby. Fragmented voids admit free energy costs and are 
thus a structural excitations. We propose that they are 
the elementary excitations enabling micro-string prop¬ 
agations, which are elementary dynamical events m- 
Also, a fragmented void is required for a micro-string 
propagation and a micro-string transports its fragmented 
void from one of its ends to the other. A few micro-strings 
can share a fragmented void producing incoherent strings 
of linear or branched geometries. Since a fragmented void 
is dominated by a smaller void with a size a fraction of 
H, in the following, we will refer to a fragmented void 
simply as a void. 

We find by inspecting hundreds of particle trajectories 
that hopping dynamics can be understood most easily 
via the motion of voids. For example, we observe that in 
an incoherent string consisting of three micro-strings at 
T = 0.32 [33] . the micro-strings propagate and reverse re¬ 
peatedly in a correlated manner most easily described by 
a single void transversing randomly along the incoherent 
string. Also, the mobile cluster of the two micro-strings 
in Fig. [7] requires two voids to properly understand. As 
collective flow subsides at low T, such descriptions be¬ 
come increasing applicable. At T < 0.32, the only exam¬ 
ples of micro-strings apparently unrelated to voids are 
loops, which are comparatively few in number and will 
be neglected for simplicity. 

We intuitively understand glassy dynamics at low T 
as dominated by string-like motions and repetitions as 
follows. A void is typically trapped by the rugged en¬ 
ergy landscape to within n -hi adjacent meta-stable sites 
with similar energies that are lower than those of the 
surrounding sites. The thermally activated discrete ran¬ 
dom walks of the void among these lower energy sites are 
equivalently repetitive string-like motions of n particles 
long. More precisely, as micro-strings can consist of two 
or more particles, the random walks of a void in general 
comprise of groups of short-range correlated hops along 
the string. After a void hops, the local energy landscape 
experienced by another void nearby may be perturbed. 
This realizes a string interaction. The motions of a col¬ 
lection of voids mutually interacting by altering the en¬ 
ergy landscape of the others generate the mobile clusters 
and thus the glassy dynamics. Recently, we have pro¬ 
posed a distinguishable particle lattice model (DPLM) 
of glass illustrating these processes in a simplified form 
m- Equation © is in particular strongly supported us¬ 
ing equilibrium data. Moreover, we have also derived a 
local configuration-tree theory of glassy dynamics based 
on such interacting voids as observed here and in DPLM 

Ba¬ 


vin. CONCLUSIONS 

In summary, string-like motions in a typical glassy fluid 
are studied by MD simulations. As T is decreased, parti¬ 
cle hops of distances comparable to the particle diameter 
are shown to dominate the dynamics. These hops and 


9 


the string-like motions they constitute become increas¬ 
ingly repetitive at low T. We propose the repetitive part 
of the motions to be /3 relaxations. At the same time, 
disruptions to the repetitive string-like motions by pair- 
interactions among strings lead to structural relaxations. 
Quantitatively, the probability P ret of hopping reversal 
reaches 73% for equilibrium samples and can be even 
higher at lower T for quenched samples. Extended an¬ 
nealing of quenched samples are found to further increase 
P re t- The rate of structural relaxation as determined 
from particle MSD is close to the rate of non-returning 
hops but is much smaller than the simple particle hop¬ 
ping rate. 

We are currently extending the present study of string 
repetition and interaction to other systems exhibiting 
string-like motions, although we expect similar results 
as these motions are remarkably similar in various sys¬ 
tems mmm- The highly repetitive nature of hopping 
motions is bound to have strong quantitative impacts on 
the overall dynamics. We believe that detailed theories of 
glassy dynamics should address or incorporate this find¬ 
ing. 

An important early model of dynamic facilitation is 
the facilitated kinetic Ising model of Fredrickson and An¬ 
dersen (FA) [9] [43]. In the two-spin facilitation version, 
a spin can only change state if two neighboring spins 


are excited as up-spins. An elementary excitation in the 
form of an up-spin has been interpreted as a fluid region 
that has a lower than average density [43] . Our identi¬ 
fication of an elementary excitation as the existence of 
a micro-string which must possess a void is consistent 
with this picture. Also, our observation of string inter¬ 
actions in which one string can be enabled or disabled 
by another string represents facilitation between excita¬ 
tions analogous to those in the FA and related kinetically 
constrained models (KCM) [IT] 

We have shown distinct T dependences of the rates R\ 
and R 2 of respectively hops and second non-returning 
hops. This can be interpreted as a difference in the ac¬ 
tivation free energies. While the displacement due to a 
single hop is around 0.9cr, a non-returning hop generates 
a larger total displacement around 1.6cr as explained in 
Sec. [V] This is consistent with an activation energy in¬ 
creasing with the displacement magnitude as proposed in 
Refs. ED [25]. 
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